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Abstract 

The time averaged and unsteady density fields close to 
the nozzle exit (0.1 < x/D < 2, x: downstream distance, 
D: jet diameter) of unheated free jets at Mach numbers 
of 0.95, 1.4 & 1.8 were measured using a molecular 
Rayleigh scattering based technique. The initial 
thickness of shear layer and its linear growth rate were 
determined from time-averaged density survey and a 
modeling process, which utilized the Crocco-Busemann 
equation to relate density profiles to velocity profiles. 
The model also corrected for the smearing effect caused 
by a relatively long probe length in the measured density 
data. The calculated shear layer thickness was further 
verified from a limited hot-wire measurement. Density 
fluctuations spectra, measured using a two- 
Photomultiplier-tube technique, were used to determine 
evolution of turbulent fluctuations in various Strouhal 
frequency bands. For this purpose spectra were obtained 
from a large number of points inside the flow; and at 
every axial station spectral data from all radial positions 
were integrated. The radialy-integrated fluctuation data 
show an exponential growth with downstream distance 
and an eventual saturation in all Strouhal frequency 
bands. The initial level of density fluctuations was 
calculated by extrapolation to nozzle exit. 

I. INTRODUCTION 

The motivation for this experiment is to provide 

benchmark data to help set-up initial conditions in the 

Computational Aeroacoustics (CAA) simulation of high- 

speed jet noise. The present work focuses on the vicinity 

of the nozzle exit. A detailed survey of the flow-field 

further downstream was presented earlier 1 . A second 
motivation for this paper is to determine the evolution of 
turbulent fluctuations in various Strouhal frequency 
bands using the present and earlier data. Such 
information can be used to validate CAA codes. 

Modem CAA methodologies, such as Direct 

Numerical Simulation 2 (DNS), Large Eddy Simulation 3 
(LES), Reynolds-averaged Navier-Stokes (RANS) 
computation 4 * * * * and semi-analytical instability wave based 
approach 5, 6 , require information on the shear layer at the 

nozzle exit. The axial and radial velocity profiles and 
intensity and azimuthal modes of turbulence fluctuations 
at the nozzle exit dictate jet development farther 


downstream. It is known that the intensity of far field 
acoustic emissions is directly related to such 
information. In fact, to obtain a reasonable comparison 
between computed and measured acoustic fields the 
initial conditions are used as adjustable parameters in 
some analysis 5, 6 . Such an arbitrary approach has been 
justified due to the lack of any experimental data. 
Furthermore, except for Direct Numerical Simulation, all 
other computation schemes selectively compute a 
portion of the turbulent fluctuations occurring within a 
range of Strouhal frequencies; thereby, imposing a low- 
pass filter function to the wide-band turbulent 
fluctuations encountered in a real flow. Validation of 
these newer computational codes requires experimental 
data not only on the time-averaged quantities but also on 
development of various Strouhal frequency fluctuations 
present in a natural jet. Once again, such experimental 
data are scarce in the available literature. Direct 
Numerical Simulation codes endeavor to compute the 
complete range of turbulent fluctuations from the low 
frequency energy containing range to the high frequency 
Kolmogorov range. At the present time such calculations 
are limited to very Low Reynolds number flow. 
Nevertheless, future maturation towards high Reynolds 
number conditions will depend on the availability of 
reliable experimental data. 

In general, the available experimental data on 
the initial conditions and development of turbulence 
fluctuations are either for low speed, incompressible 
jets 7, 8 ' 9 ‘ 10 or for low Reynolds number high speed jets 11, 
12 . There is a void of information when it comes to 
practical, high Reynolds number, high-speed jets. The 
problem lies with the traditional experimental 
techniques. Until now most unsteady turbulence 
measurement were performed using a hot-wire probe 
which is suitable for low speed flows where density and 
temperature variations are insignificant. In spite of the 
ambiguity in signal analysis, hot-wires are also used in 
low Reynolds number supersonic flows where 
aerodynamic load on the wire is small. In high Reynolds 
number flows hot wires are prone to breakage, which 
ultimately limits their application. The particle based 
optical techniques, such as, Laser Doppler Velocimetry 
(LDV) and Particle Image Velocimetry (PIV) are 
primarily suitable in providing time-averaged 
measurements. The present work relies on a Rayleigh 
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scattering based point measurement technique that has 
been demonstrated to provide density fluctuation 
spectrum in the compressible flow regime 1 . The 
technique is under development to provide velocity 
fluctuation spectrum 13 , which will provide a substitute 
for hot-wires at high-speed flows. The present work is 
based on the proven ability of the Rayleigh scattering 
technique to provide time averaged and unsteady air 
density measurements. Basically, a laser beam is passed 
through the flow and light scattered by gas molecules are 
collected and analyzed to simultaneously obtain 
information on velocity, temperature and density. Air 
density, the quantity of interest, is proportional to the 
intensity of scattered light. 

II. EXPERIMENTAL SETUP 

The experiments were performed in a small jet facility 
at the NASA Glenn Research Center. A convergent 
nozzle was used to produce a Mach 0.95 plume and 
convergent-divergent nozzles were used for the Mach 
1.4 and 1.8 plumes. All three nozzles had a nominal 
equivalent diameter of D = 2.54 cm. 

Rayleigh scattering principles and a version of 
the measurement technique have been described in detail 
in references 1 & 14. The current optical arrangement 
around the jet facility, shown in Fig. 1, is somewhat 
different. It provides improvements of increased laser 
power and accessibility to the flow region near the 
nozzle exit. The major difference between the present 
and the earlier experimental setup is the optical 
arrangement near the nozzle. The earlier setup used a 
larger diameter beam that was focused at the probe 
volume and passed at 45° to the jet axis. Interception of 
a part of the beam by the nozzle rig limited the closest 
measurement point to be about 1.2D downstream of the 
exit. The present setup, on the other hand, uses a narrow 
beam perpendicular to the nozzle axis. This arrangement 
allows measurements as close as 0. ID from the exit. In 
brief, a polarized, narrow diameter laser beam from a 
continuous wave, frequency doubled NdiVCU laser 
(532nm wavelength) was passed perpendicular to the jet 
axis. The light scattered by the air molecules were 
collected vertically below, at 90° from the incident light 
direction. The f/# 3.3 collection lenses focused the 
scattered light on the face of a 0.55mm diameter 
receiving optical fiber. The fiber diameter and the 
magnification ratio of the collection optics fixed the 
probe volume length at 1 .03mm. The beam waist was 
about 0.1 6mm 2 in cross-section. In effect, light scattered 
from this small length of the beam was collected by the 
receiving fiber. The collected light was then split into 
two parts and measured using a pair of photo multiplier 
tubes (PMT) and photon counting electronics (not shown 
in Fig. 1). This opto-electronics part is identical to that 
of reference 1. Both the laser source and the receiving 
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optics were placed on a X-Y traverse that allowed the 
probe volume to be moved from point to point in the jet 
plume. An important concern in Rayleigh scattering 
setup is dust removal from the air streams. In the present 
setup the primary air was cleaned using air filters that 
blocked particles above a micron diameter. In addition, 
an 8" diameter low speed 20m/s) co-flow around the 
primary jet was maintained using a second filtered air 
source to avoid particles through the entrained air. 

Fundamentally, for a constant molecular 
composition air and a fixed optical system, the intensity 
of scattered light is directly proportional to the local air 
density p. Since light intensity is measured in terms of 
the number of photoelectrons N counted over a given 
time interval At, the following relationship holds 

N-(ap + b)At (1) 

Here a and b are constants determined through an in situ 
calibration. The calibration process was the first step in 
density measurement. It was performed in the plume of the 
baseline circular convergent nozzle operated in the Mach 
number range of 0 to 0.99. At each operating condition the 
photon arrival rate was counted over one second duration 
and the jet density is calculated using isentropic relations. 
Subsequently, a straight line was fitted through the data to 
determine the proportionality constants a and b. Since two 
PMT and photoelectron counters were used, two sets of 
calibration constants ai, bi and a 2 , b 2 were calculated. The 
photoelectron counting was performed over a series of 
contiguous time bins. The statistical mean, when 
multiplied by the calibration constants, provided time- 
averaged air density. The time-averaged data are quite 
accurate: absolute density numbers are found to be 
repeatable within ±1% of quoted values. 

The 2 PMT technique was necessary to measure 
fluctuation spectrum. A simpler approach would be to 
use a single PMT, perform a single series of photon 
counting, and to take a Fourier transform of the series. 
That method, however, is affected by electronic shot 
noise inherent in any optical measurements. To reduce 
the effect of shot noise, the collected Rayleigh scattered 
light was split into two parts, measured using two PMTs 
and photon counting electronics and finally the two 
series of photon counts were cross-correlated. Since, 
shot noise emitted by individual PMTs is uncorrelated, 
contribution from this source is significantly reduced. 
The simultaneous photoelectron counting produced two 
series of data Nh and N 2 , (i = 0, 1,2,.... n-1). The average 
values from each of the time series were subtracted: N'i,= 
Nh - Ni av , N' 2i = N 2i - N 2av and a cross spectral density 
P was calculated from individual Fourier 

transform, F , (/), F N , (/) calculated at discrete frequency 

bands f ) : 




Superscript * in the above equation indicates complex 
conjugate. The density fluctuation spectra is calculated 
using appropriate calibration constants a, and a 2 for the 
two photomultiplier tubes: 



f) 

a l a 2 {\t) 2 


( 3 ) 


Usually two long records, typically of 524,288 data points, 
were collected from multiple segments of 16,384 data 
strings. The latter is the maximum number of contiguous 
counts delivered by the photon counters. The Welch 
method of modified Periodograms 15 was used to calculate 
the cross-spectral density. Each long record was divided 
into small segments of 512 data points. The adjacent 
segments were overlapped by 50%. The modified 
periodograms of corresponding segments from the two 
PMTs were calculated and then used to determine local 


estimates of cross-spectral density. All local estimates 
were averaged to obtain the final cross-spectral density. In 
all spectral data presented in this paper the Nyquist range 
was divided into 256 frequency bins. 

For validating the shear layer thickness data, the 
time-averaged axial velocity profile in the Mach 0.95 jet 
was estimated using a single hot-wire probe. The 
combination of a thin shear layer close to the nozzle exit 
and a relatively long laser probe length introduced 
additional modeling complications which needed to be 
validated by an independent means. The hot-wire 
technique, while providing adequate spatial resolution, 
suffers from many logistical as well as signal analysis 
problems. Sensor survivability is by far the most serious 
logistical problem. In compressible flow analysis of hot- 
wire signal is complicated by its dependence on air 
density, temperature as well as velocity. While accurate 
measurements are possible by using multiple overheat 
ratios and multiple sensors, they are not straightforward 
and may encounter other difficulties. (For example, use 
of multiple sensors would again run into the spatial 
resolution problem). At transonic speeds the dependence 
of hot-wire signal on flow Mach number brings 
additional complications. In this experiment a single hot- 
wire was simply operated with a single overheat ratio. 
The probe was placed at the exit of the nozzle, and the 
voltage output was calibrated against the jet velocity. 
This calibration was used to measure the velocity in the 
shear layer. For an estimate of the shear layer width this 
procedure was deemed adequate. No attempt was made 


to perform such measurements at the supersonic 
conditions. 

The hot-wire probe was inserted at an angle 
from the low-speed side so that interference with the 
flow was as small as possible. The profiles were 
obtained in the range, 0.1 < x/D < 1, with a probe- 
traversing unit under automated computer control. 
Several attempts to measure the ‘initial profile’ just 
downstream of the nozzle (x/D « 0.03) resulted in probe 
breakage. The sensor typically broke when reaching the 
low speed edge of the shear layer. The reason for this 
has remained unclear. However, the profiles at this 
location could be obtained up to a jet Mach number of 
0.7. Some of these results will be discussed in the text. 

III. RESULTS AND DISCUSSION 
A. Modeling of time averaged density data: Shear 
layer parameters are traditionally expressed through 
axial velocity profiles. Since, the present work 
represents jet flow in terms of air density variation, a 
need for suitable modeling and the use of Crocco- 
Busemann’s equation to convert the modeled velocity 
profile to density profile became necessary. An iterative 
process was setup where the first step was to assume a 
velocity profile (equation 4) with a single variable b that 
determined shear layer thickness. For an estimated value 
of b the density profile was calculated using Crocco- 
Busemann equation (equation 6) and compared with 
experimental data. The estimate of the b parameter was 
improved for the next iteration until the difference 
between the modeled and measured profiles was 
minimized. For radial profiles obtained close to the 
nozzle exit, the averaging effect caused by a relatively 
long probe volume had to be included in the model. 
Instead of directly using Crocco-Busemann (equation. 6) 
an integrated form (equation 8) had to be used. The 
details of the modeling process are discussed in the 
following. 

For velocity profile the patched half-Gaussian 
profile was chosen, from a large number of alternatives 
(Michalke 16 ), based on its wide application in the 
existing literature 5 ’ 12, 17 ' 18 . 

U j Jr-h ) 2 1 

— = expj-m2 > for r>h 

Uj { b )\ J (4) 

= 1 for r <h 

Where, r is the radial distance, h is the height of the 
potential core at the measurement station, U is time- 
averaged velocity, Uj is the centerline velocity and b is 
the half width of the mixing layer (Fig. 2)**. The above 

** In incompressible flow it is customary to describe the 
boundary layer at the nozzle exit through momentum 
thickness: 
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expression is for top-hat profiles and is valid till the end 
of the potential core: well within the range our interest. 
The potential core height h is not an independent 
parameter, as it can be determined from momentum 

conservation, once b is known. 

00 

| pU 2 rdr = const, at all axial stations (5) 

o 

The normalized density profile is obtained from the use 
of Crocco-Busemann law 19 : 


-c. 

M 

2 

„ u „ 


+ C> b Cj 

u i 


u, 


\ J J 

J 


y—l 9 T ( y-1 f\ 

C.= r - — M 2 , C 2 =^- 1 + 4- — M 2 , (6) 

2 T\ 2 J 

c 3 — Q — C 2 + 1 ■ 

Here, M is jet mach number, T a ambient temperature and 
T r plenum temperature. The validity of the Crocco- 
Busemann law in jet plume has been checked 
computationally in the past 6 . In summary, for given 
Mach number, plenum and ambient temperatures, the 
density profile at a given axial station is uniquely 
dependent on the parameter b, the mixing layer half 
width. Figure 3 presents some sample radial density 
profiles and best- fit model profiles. Usually, the plenum 
and ambient temperatures, T a and T r were adjusted 
within ±2°K at the beginning of the iteration to match 
the measured ambient and core density. These two 
temperatures were found to drift slightly over the 
duration of the experiment. 

Modeling of the experimental density profile 
required an additional consideration for the relatively 
large probe volume length (28 = 1.03mm or slightly 
more than 4% of jet diameter). The measured density 
data is an average of the distribution over the probe 
volume length, 28: 




dr , 


which, for the above half-Gaussian profile can be shown 
to be 0 inc = 0.31 1775b. For compressible flow the 
definition of momentum thickness includes density 1 : 



Jrf i 

Pi tu 



dr , 


and thereby making it dependent on the plenum and 
ambient density. The compressible version is infrequent 
in the literature and, therefore, the mixing layer is 
described only through the half width parameter, b in the 
present paper. 


( 

Pj 

\ r+S 

= — f 

\ 

p J 


28 J 

' meas r o\ 

,P) 


The averaging effect is insignificant for axial positions 
x/D > 1 , where shear layer is much thicker than the 
probe length. However, the thin shear layer close to the 
nozzle exit warrants special consideration. Even when 
the shear layer thickness is smaller than the probe length, 
radial traversing in a fine increment (minimum of 
0.025mm in the present experiment) produces a smeared 
profile as described in the above equation. Figure 4 
presents a schematic of the smearing process. Since, the 
probe volume length is known with high degree of 
fidelity, the true steeper profile was recovered from the 
smeared data. (As mentioned earlier, the probe volume 
length was tightly fixed by the magnification ratio of the 
receiving optics and the aperture imposed by the 
receiving optical fiber core diameter.) In the modeling 
approach, the probe-length-averaged model profile was 
calculated for various values of the shear layer half- 
thickness b and checked for best fit with the 
experimental data. The patched profile of equation 4 
required that the integration of equation 7 be made over 
three separate regions: 



+ (C 2 - 1 ){r- 5-lt) ] + C 2 for h + S> r >h - 5 
= 1 for r <h—S 

, r . — — ( r — S — h ) r — —(r + 8—h) 

where, a = Vln2 , f - Vln2 , 

b b 

Ci b C,b 

4 5 VM2 

For measurement stations close to the nozzle exit the 


iterative scheme involved numerically computing 
equation 8 for estimated values of b and comparing with 
experimental data. In Figure 3 solid lines show best- fit 
profiles. Once the best- fit value of b parameter is 
established the corrected density profile, without 
smearing effect, is back calculated from equation 6. 
Figure 3 shows the corrected profiles by chain lines. By 
comparing the two lines it can be seen that the averaging 
effect is substantial for profiles at x/d = 0.1 and minimal 
at x/D = 1.0. Usually the modeled profile shows 
excellent fit to the measured density data. The optimum 
value for the shear layer half width, b can be estimated 
with fairly high accuracy ±2% for measurement stations 
x/D > 1. The uncertainty increases to ±10% for the 
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station closest to the nozzle exit. Additional problems 
arise in the supersonic plumes, such as Fig. 3(d), where 
weak shocks cause significant deviation from the 
modeled top-hat profile. Nevertheless, weak shocks 
primarily affect the potential core, leaving the shear 
layer nearly intact. 

Figure 5 presents the variation in shear-layer 
half-width, b, with downstream distance measured from 
a large number of radial traverses and modeling 
procedure described above. There are several pieces of 
information available from this figure. First, the shear 
layer grows linearly with downstream distance for all 
Mach number conditions. This information can be used 
to determine the radial velocity V distribution, via the 
use of the axisymmetric continuity equation: 

d(pv) | pv | d(pu ) _ 0 

dr r dx 

Second the growth rate progressively reduces with 
increasing Mach number. Both of these follow earlier 
observations 12 ' 16 . Third, from the linear least square fit, 
the shear layer half-width at the nozzle exit is estimated 
as 0.0013, 0.0039 & 0.0144 for Mach 0.95, 1.4 and 1.8 
jets respectively. The increase in shear layer thickness in 
the supersonic plumes may be due to interaction between 
nozzle boundary layer and weak shocks inevitably 
present in the diverging part of nozzles. 

B. Validation from hot-wire measurement: An 

indirect route has been used in the present work to 
determine the shear layer parameters through density 
measurement and a model fitting procedure. Any other 
technique for accurate measurement of these parameters 
is presently unavailable. The problems associated with 
the hot-wire technique have been discussed in section II. 
Two velocity profiles obtained by the hot-wire technique 
along with the modeled profiles are shown in Fig. 6. The 
modeled profile is that of equation 4 with half thickness 
b same as that calculated from density data. The 
agreement between the hot-wire profile and the model 
profile is excellent in Fig. 6(a) and reasonable in 6(b), 
which provides some confidence in the indirect route 
used to determine the shear-layer thickness, b. 

The hot-wire was also used to measure shear- 
layer profiles extremely close to the nozzle exit, x/D = 
0.03, for low Mach number subsonic jets. These 
measurements provide a glimpse on the initial boundary 
layer state that may be critical in computational efforts. 
Figure 7 presents the incompressible momentum 
thickness 0 inc (see footnote) calculated directly from the 
measured velocity profiles. As stated in § II, these data 
could be obtained up to a jet Mach number of about 0.7. 
The data point for Mach 0.95 jet was obtained from 
Rayleigh measurement of shear layer half width b and its 
relation to the momentum thickness, 0i nc = 0.3 1775b. 
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The momentum thickness was found to vary inversely 
with the square root of the Reynolds number, based on 
exit diameter, Re D . The chain line in Figure 7 is an 
empirical fit: 0i nc /D = C/Re D 12 . The value of C that is 
found to provide a good fit is 0.77. From shape factor 
values and trends of the data the boundary layer is 
inferred to be ‘nominally laminar’. Note that the shear 
layer is very thin. The momentum thickness is only 
about 0.1 percent of D (D/0 « 1000). Additionally, that 
the Rayleigh measurement of momentum thickness 
followed the same trend observed from hot-wire data 
provides confidence in the modeling process. 

C. Density fluctuations spectra and evolution of 
instability waves: The density fluctuation spectra are 
measured using a two-PMT technique described earlier. 
Figure 8 presents spectra measured at various points 
along the lip line of a Mach 1.8 jet. Note that the 
frequency f is non-dimensionalized to Strouhal number 
Sr = fD/Uj and spectral energy in the ordinate is non- 
dimensionalized by (pj - p a ) 2 . The spectral data have a 
fundamental noise floor and random uncertainty. The 
former is measured by turning the primary jet flow off, 
that is, under no density fluctuation condition and the 
latter can be seen as superimposed randomness on each 
spectral plot. The mean square of this randomness is 
about 5x1 O' 6 (kg/m 3 ) 2 . This randomness dominates spectra 
from axial stations with very small density fluctuation, 
such as the first measurement station at x/D = 0.1. Further 
downstream, not only does the spectral energy increase, 
but also the spectral peak shifts from Sr = 2.5 to 0. 1 . The 
straight line in Figure 8 shows the existence of the inertial 
range in turbulent fluctuations with -5/3 slope. In the shear 
layer of a turbulent jet, high frequency fluctuations are 
generated by two means. First, high frequency Kelvin- 
Flelmholtz instability waves are generated in the thin shear 
layer adjacent to the nozzle exit. Second, energy cascading 
extracts energy from the low frequency instability waves 
and creates the inertial sub-range. This second means is 
established further downstream from nozzle exit. Figure 8 
confirm s the emergence of this scale around x/D = 7 in 
Mach 1.8 jet. 

To determine evolution of instability waves a 
large number of density fluctuation spectra were measured 
close to the nozzle exit. Additional data, taken further 
downstream of the nozzle exit and available from previous 
effort 1 , were also utilized. In the present experiment a 6X6 
grid covering typically 0.1 < x/D < 2.0 in the axial 
direction and the shear layer width in the radial direction 
was used. The earlier data typically used a 10X10 grid that 
covered 1.2 < x/D < 12.0 and shear layer width. The exact 
number of stations was somewhat different for different 
Mach number jets. At a given axial position band-passed 
density fluctuations were integrated for all radial locations 
similar to that presented by Gaster, Kit & Wygnanskii 20 for 



a low speed mixing layer. Following the notation of 
equation 3 the process can be described as: 


r ou, 



where r in is the radial position for the innermost 
measurement station and r out is the outermost. Both r in and 
r out were varied to cover the shear layer thickness. For 
example at x/D = 0.1 r in = 0.45 and r out = 0.53 while at the 
end of the potential core r in = 0 and r out = 1 ,0D. The band- 
passed fluctuation at a given Strouhal frequency was 
determined over a 195.3 FIz. wide band. The frequency 
width was fixed from 256 frequency bins covering 50 kFlz 
range. The latter in turn was fixed from the time duration 
(At = 1 Omicro-seconds) for the individual photon counting 
bin. Figure 9 presents a comparative study of the evolution 
of the band-passed turbulent fluctuations in the three Mach 
number jets. Note that the band-passed fluctuations were 

normalized as p' b /( p • - p a ) for a uniform comparison in 

Figure 9. Similar turbulent evolution data for low 
Reynolds number supersonic jets were presented by 
Morrison & McLaughlin 12 and for low speed mixing 
layers by Gaster, Kit and Wygnanskii 20 . 

Evolution data of Fig. 9 is best described through 
a comparison with the normal mode representation 16 of 
hydrodynamic instability waves. Note that stability 
analysis for the present flow conditions has not been 
performed; only a qualitative comparison with similar, 
past analyses has been made in the following. Fluctuations 
can be expressed through any scalar (density, temperature, 
pressure) or vector (velocity) quantity in the normal mode 
representation. For example, the instantaneous density 
p (x, r, <)), t) is expressed as a sum of fluctuations occurring 
at frequency co, azimuthal mode m, and axial wave number 
a: 

p' (x,r,<f>,t) = 


IE p(r, x, co) exp j / J adx + m (j)-cot 


where p (x, r, co) is the eigenfunction and a is eigenvalue 


of the problem. The latter is a complex number, a = a r - 
ia; whose real part is related to the phase speed and the 
imaginary part dictates local growth or decay rates of the 
instability waves. The radian frequency co is related to the 
Strouhal frequency Sr as: 


2/rf7 Sr 

co- 

D 


( 12 ) 


The integration in the wave number a accounts for the 
cumulative change in the growth rate of the instability 
waves due to the divergence of the shear layer 7 . The 


unsteady data presented in this paper can be related to 
above analytical expression. The spectral data resolves 
fluctuations in various Strouhal frequencies. Flowever, 
resolution in the azimuthal modes m was not possible. The 
spectral data, therefore, is a summation over all azimuthal 
modes. The expression for the density spectra in equation 
3 is effectively the following: 

co r Gv ' 

JP^7(x,r,Sr)= Y, L o(x, r, co) expt i\ J a dx + m (j) - cot 


- p' {x,r,co)QX^\i | 


a dx - cot 


The radialy integrated band-passed fluctuation (equation 
10) approximates the following: 

00 f f x 

p' bp (x, Sr) -2n\ p ' (x, r, Sr) expj i \adx-cot > rdr 
o l Vo J } 


i | a dx 


2nU Sr 


Multiple observations can be made by relating 
the above equation to the evolution plot of Figure 9. 
Since the evolution data are plotted in a semi-log scale a 
linear increase with the axial distance reflects an 
exponential growth 12 as expected from equation 14. This 
exponential growth of fluctuations is seen in all 
frequency bands close to the nozzle exit. From the 
hydrodynamic stability analysis it is known that the slow 
divergence of the jet shear layer progressively reduces 
the amplification rate, thereby saturating (and decaying 
further downstream) fluctuation amplitudes 5 ’ 6 ’ 1 . The 
experimental data follows this trend. Additional 
observations that follow expectations from stability 
analysis are the lower the Strouhal frequency the higher 
the peak amplitude, and the axial location of the peak 
occur further downstream. An examination of Sr = 0. 12 
data in Mach 1.8 jet shows that the fluctuation level 
increases from a value of 0.0005 to 0.05: a factor of 100! 

An estimate of the fluctuation levels at the 
nozzle exit was obtained by linearly extrapolating the 
evolution curves. The values are presented in Table I. 
Since the extrapolated values do not show a consistent 
trend with increasing Strouhal frequency, a range is 
specified. The range also reflects the random uncertainty 
associated with the calculation and extrapolation 
process. The initial fluctuations levels are highest for the 
lowest Mach 0.95 jet and progressively decreases with 
increasing Mach number. 

The uncertainty in the evolution data primarily 
arises from three sources. First, the random uncertainty 
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in density spectra measurement: estimated to be about 
5x1 0‘ 6 (kg/m 3 ) 2 . However, the summing process over all 
radial stations is expected to reduce this uncertainty. 
Second, the noise floor in the spectral data, estimated 1 to 
be 4x1 O' 6 (kg/m 3 ) 2 , and provides a bias error in the 
evolution plot. The summing process increases the bias 
error. For example, sum over 10 spectral estimates 
makes the bias error 4xl0" 5 (kg/nr 3 ) 2 . The third source is 
the limited radial extent (r out < r < r in ) over which density 
fluctuations are measured. Ideally, fluctuations from the 
convective eddies extend to infinity, i.e., r out = cc; 
however, the fluctuations amplitudes rapidly decay 
beyond the turbulent flow boundary 6 . The uncertainty 
due to limited radial extent is difficult to establish, yet it 
is expected to be lower close to the nozzle exit and 
higher farther downstream. It is estimated that the 
uncertainty in the band-passed data is in a range of 5% 
to 20% of the presented values. The highest uncertainty 
appears close to the nozzle exit and in Mach 0.95 plume, 
while the lowest uncertainty is for measurement stations 
farther away from nozzle exit and in Mach 1.8 plume. 
Since the difference in air density between the jet plume 
and ambient air increases with increasing Mach number 
the spectral estimates become more reliable at higher 
Mach number conditions. 

An important parameter, that the present 
experiment could not provide, is the distribution of 
turbulence fluctuations in various azimuthal modes. 

Such data would have required multiple measurements 
in various circumferential positions, which is 
unavailable at the present time. Nevertheless, DNS 
calculations 2 , and earlier near field microphone 
measurements 10 indicate that the first few azimuthal 
modes (up to ±2) account for bulk of the turbulent 
fluctuations. 

IV. CONCLUDING REMARKS 

The initial condition data presented in this paper, along 
with the noise source identification in reference 1 and 
the time-average==d velocity and temperature data 
presented in reference 2 1 is expected to provide 
sufficient information for validation of various 
Computational Fluid Dynamics (CFD) and 
Computational Aeroacoustics (C AA) codes. All of these 
measurements were performed in the same plumes of 
Mach 0.95, 1.4 and 1.8 jets. The initial status of the 
shear layer at nozzle exit is surmised in table 1. It is 
hoped that the advancement of the Rayleigh scattering 
technique will enable further insight into the turbulent 
character of the jet including evolution of velocity 
fluctuations and 2-point space-time correlations. 
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Table I. Summary of nozzle exit conditions: 

Specific heat ratio, y = 1.4; Nominal ambient and total temp., T 0 = 300°K 


Plume Mach number 

0.95 

1.4 

1.8 

Reynolds number Ren 

0.66xl0 6 

1.16x10 s 

1.88xl0 6 

Shear layer half width at nozzle exit, b/D = 

.0031 

.0039 

.0144 

Rate of growth, db/dx = 

0.1104 

0.0873 

0.0654 

Initial fluctuation level over 1 95 Hz. band 
P bp/(Pi - Pa) 

0.001 -0.0014 

.0006 - 0.0007 

0.0004 - 0.0006 



Fig. 1 . Photograph of the facility 


NAS A/TM— 2004-2 1 2392 




Distance from centerline, r/D 

Fig. 3. Time averaged density profiles from indicated axial stations and Jet Mach numbers. + +: experimental data; 
: probe-volume-averaged profile (equation 8); : corrected profile without averaging effect. 



square linear fit (described in legend). 
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Fig. 8. Density fluctuation spectra from Mach 1.8 jet 
along r/D = 0.48 and at indicated axial positions. 
Successive plots are shifted by a multiplication factor of 
2 . 
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